C/O white dwarfs of very low mass: 0.33-0.5 M 
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Abstract. The standard lower limit for the mass of white dwarfs (WDs) with a C/O core is 
roughly 0.5 Mq. In the present work we investigated the possibility to form C/O WDs with 
mass as low as 0.33 Mq. Both the pre- WD and the cooling evolution of such nonstandard 
models will be described. 



1. Introduction 

It is by now a well-known and commonly accepted result of the stellar evolution theory that stars 
of low and intermediate mass will end their lives as white dwarfs (WDs). Since the pioneering 
papers by Chandrasekhar (1931) and Stoner (1932), it has been known that there must be an 
upper limit of the mass of WDs. The current value for such a limit, now universally known as 
Chandrasekhar's mass, is M^= 1.456 (2/fi e ) 2 M (where \i e is the electronic mean molecular 
weight, Hansen, Kawaler & Trimble 2004). Below this critical mass there are essentially three 
different families of WDs concerning the chemical composition of their cores: the O/Ne/Mg, C/O 
and He. The upper mass limit for C/O WDs, which is still subject to some debate, since it is the 
aftermath of the evolution of a star with initial mass slightly lower than M up , the lowest mass 
for star that succeeds in igniting carbon burning before the asymptotic giant branch (AGB), it 
is around 1.05 MQ(Dominguez et al. 1999, Weideman 2000, Prada Moroni & Straniero 2007, 
Catalan et al. 2008, Meng, Cheng & Han 2008). 

On the other hand, the value for the minimum mass of WDs with a C/O core is commonly 
believed to be around 0.5 M (Weideman 2000, Meng et al. 2008). The theory of stellar 
evolution predicts that an isolated single star with mass lower than about 0.5 Mq, the exact value 
depending on the chemical composition (0.46 and 0.50 for Z=0.04 and Z=0.0001, respectively, 
Dominguez et al. 1999, Girardi 1999), do not ignite the helium-burning and die as WDs with an 
helium rich core. For these stars, during the red giant phase the electronic thermal conduction 
and the emission of neutrinos, which cool the core, prevail on the heating due to the contraction 
caused by the accretion of fresh helium processed by the hydrogen-burning shell. This is the 
main reason why WDs of mass less than about 0.5 M are commonly believed to have a He-core. 

On the other hand, as early as 1985 Iben &; Tutukov, in a very beautiful and instructive 
paper, showed that a star with an initial mass of 3 Mq evolving in a close binary can produce 
a remnant of 0.4 Mq with a sizeable C/O core. A similar result was obtained 15 years later 
by Han, Tout &; Eggleton (2000), who computed the evolution of close binary systems with 
several different parameters. They showed that in a binary system with the primary star of 



mass Mi= 2.51 Mq, the mass ratio q=Mi/M2=2 and the initial period P=2.559 d, the primary 
star succeeded to ignite the helium-burning. They did not followed the evolution of the remnant 
C/O WD, since their code met numerical problems when the primary star was as low as 0.33 
M with a C/O core of 0.11 M . 

In this paper we will investigate both a possible evolutionary scenario able to produce C/O 
WDs of low mass, that is M < 0.5 Mq, and their cooling evolution. We will show the results of 
detailed evolutionary computations performed with a full Henyey code able to follow consistently 
the evolution of stars from the initial pre-main sequence to the final cooling phase of WDs. 

2. The red giant phase transition 

The stellar models showed in the present work have been computed with an updated version 
of FRANEC (Prada Moroni & Straniero 2002, Degl'Innocenti et al. 2008), a full Henyey 
evolutionary code. We adopted a metallicity Z=0.04 and an initial helium abundance Y=0.32 
suitable for the very metal-rich stars belonging to some open galactic cluster, such as NGC6791. 
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Figure 1. Mass of the He-core at the tip 
of the RGB (diamond and dashed line) 
and mass of the hydrogen exhausted core 
at the first thermal pulse on the AGB 
(circles and solid line) as a function of 
the initial mass for stars with Z=0.04 
Y=0.32. 



Figure 2. Mass of the hydrogen 
exhausted core as a function of the age 
for a star ofM=2.3M , Z=0.04 Y=0.32. 
The circles mark: the central hydrogen 
exhaustion, the base and the tip of the 
RGB and the central helium exhaustion. 



The figure Q] shows the mass of the hydrogen exhausted core M# at the first thermal pulse 
on the AGB (upper curve) and the mass of the helium core M# e at the tip of the red giant 
branch (RGB tip, lower curve) as a function of the initial mass. As it is well known, M# at 
the first thermal pulse is nearly constant around 0.55 Mq, the exact value depending on the 
chemical composition, for initial masses lower than 3 Mq. On the other hand, the behavior 
of the M.He a t the tip of RGB as a function of the initial mass is much less smooth than 
the previous one, in particular a deep and sharp minimum is present around 2.3 Mq for this 
chemical composition (Castellani, Chieffi, Straniero 1992). Such a behavior is the consequence 
of the physical conditions present in the helium-core at the onset of the 3a nuclear reaction. 

For initial mass lower than 1.5-1.7 Mq, an electron-degenerate core develops at the beginning 
of the red giant phase, and the onset of the 3a occurs in a nuclear runaway, the so-called helium- 
flash (Hoyle & Schwarzschild 1955). For these stars, the M# e at the tip of the RGB is almost 



constant, about 0.46 M for metal-rich stars and 0.50 Mq for metal-poor ones, and is essentially 
produced during the hydrogen-shell burning phase. 

Prom this point on, the larger the initial mass, the weaker the degeneracy of the core in the 
red giant phase, thus the weaker the He-flash and the lower the mass of the helium-core required 
for the 3a burning. When the initial mass is as high as about 2.3 Mq, the helium-burning ignites 
quiescently, without a flash. This explains the steep decrease of the M# e at the tip of the RGB 
between 1.8 and 2.3 Mq. For higher masses, the core in the red giant phase is not degenerate 
any more. From this point on, the larger the initial mass, the larger the M.jj e a t the tip of the 
RGB, as the larger the convective core during the previous main sequence phase. Thus, the 
minimum mass of the helium-core M# e at the onset of the 3a occurs in the transition between 
these two regimes, around 2.3 M , reaching a value of 0.315 Mq. 

Notice that, as previously stated, this deep minimum is almost canceled at the first thermal 
pulse. In fact, the hydrogen burning shell continues to process hydrogen and move outward, 
then increasing M#, during the central helium-burning phase. Moreover, it is a well established 
result of the theory of stellar evolution that the lower the mass of the He-core at the onset of 
the 3a, the fainter the star and the longer the central helium-burning phase. Thus, as the initial 
mass of the helium-rich core decreases the hydrogen-burning shell will work for a longer time 
resulting in an larger AM//. This is the reason why star with initial mass lower than 3 M enter 
the AGB phase with nearly the same M# (Dominguez et al. 1999). 

The first paper devoted to the study of the evolutionary characteristics of stars belonging to 
this sharp and deep minimum in the mass of the He-core at the tip of the red giant phase was 
Iben (1967). Since then, many studies have been focused on this transition, named by Renzini 
& Buzzoni (1983) red giant phase transition (hereafter RGB phase transition), but a detailed 
discussion of such an important feature is well beyond the aim of the present paper. We will 
dwell only upon a feature of the evolution of a star with M= 2.3 Mq important to our purpose. 
Figure [2] shows the evolution of the mass of the hydrogen-exhausted core M# for this transition 
mass. The circles mark the following evolutionary phases: the central-hydrogen exhaustion, the 
base and the tip of the RGB, respectively, and finally the central-helium exhaustion. As one can 
easily see, when the star exhausts its central hydrogen, the mass of the He-core is already 0.225 
Mq, which is more than the 70% of the value at the RGB tip. In fact, a sizeable convective core 
was developed during the preavious main sequence phase. When the star reaches the base of 
the RGB, the mass of the helium core is 0.254 Mq, which is about the 80% of the mass required 
to onset He-burning. This means that, since a large fraction of the helium-core required for the 
ignition of the He-burning is already developed at the beginning of the red giant evolution, a 
star in the RGB phase transition, that is around 2.3 M , that undergoes a very strong mass 
loss during the red giant phase can produce a C/O WD with a mass significative lower than 0.5 

Mq. 

3. C/O white dwarfs of very low mass 

In order to check such an idea and provide models of C/O WDs with mass lower than 0.5 M , 
we computed the evolution at constant mass of a star with M=2.3 Mq until the red giant phase. 
Then at about logL/Lo=1.34, when the He-core was M# e = 0.2569 M , we turned on the mass 
loss, with a constant rate of the order of 10~ 7 Mq yr _1 . We stopped the mass loss once the 
desired final mass was obtained and then we followed the next evolution until the final cooling 
phase. 

The figures [3] - [7] show the evolutionary tracks in the HR-diagram for different values of the 
final mass, from 0.48 Mq to 0.33 Mq, from the red giant to the final WD phase. The model 
showed in figure EJ with a total mass of 0.48 Mq, succeded to ignite He-burning and become 
a WD with a C/O core whose mass fraction is about 89 %, where we defined the core edge to 
be the point where the helium abundance became larger than 50%. A very similar evolution is 
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Figure 3. Evolutionary track of a star 
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followed by the model with 0.461 Mq, as can be seen in figure HI In this case, the resulting WD 
has a C/O core of 85% of the total mass. 




Figure 5. Same as in figure [3] but for a Figure 6. Same as in figure [3] but for a 
remnant of mass M=0.43 M . remnant of mass M=0.38 M . 

The evolution becomes more complex and difficult to compute for lower masses. The figure 
[5] shows the evolution of the model with 0.43 M . The star onsets the 3a nuclear reaction and 
after the central He-burning phase it experiences a series of He-shell thermal pulses: during 
each flash the star describes a loop in the HR diagram. After the end of the thermal pulses 
the helium-burning shell continues quiescently and the model moves toward higher effective 
tempearures. Then the star, approaching the cooling track, experiences a few strong hydrogen- 
shell flashes, more specifically CNO-flashes. During these episodes, the hydrogen rich outermost 
layer is progressively eroded, until the star can cool down as a WD with a C/O core of 79% 
of the total mass. The evolution of the model with 0.38 M , shown in figure El presents again 
both the He-thermal pulses and the CNO-flashes. In addition, after the end of themal pulses, 
when the model moves toward the blue it experiences a late thermal pulse. The resulting WD 
has a C/O core of 50% of the total mass. 

Notice that the lower the mass of the remnant, the lower the mass fraction of the C/O core, 
thus the larger the mass of the helium-rich buffer. In standard C/O WD, the thickness of the 
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Figure 7. Same as in figure [3] but for 
a remnant of mass M=0.33 Mq. The 
model does not succeed to ignite 3a and 
eventually cools as a He-core WD. 
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Figure 8. Same as in figure [3] but for 
a remnant of mass M=0.33 Mq. In this 
case the model ignites He-burning and 
eventually cools as a C/O-core WD. 



He-rich layer is of the order of a few percent (1-2%), while in the present models of very low 
mass C/O WDs it can reach the 50% of the total mass. 

Notice that the helium-shell thermal pulses we showed above are not the same that occur in 
AGB stars. In fact, in that case the thermal instability is the consequence of the accumulation 
of a critical mass of helium accreted by the quiescent hydrogen-burning shell. While in this 
case, they are self-excited relaxation oscillations, where the compression is due to the relaxation 
of the star after a previous contraction. At variance with thermal pulses in AGB, these pulses 
affect the entire star, with large oscillations in central temperature and density, as early shown 
in the pioneering paper by Iben et al. (1986). 

The evolution of the remntant with mass M= 0.33 Mq deserves a detailed discussion. In 
fact, adopting the same mass loss rate along the red giant phase of the progenitor star of 2.3 
M as in the previous cases, the 0.33 M does not experience 3a burning and becomes a WD 
with a helium core. Figure [7] shows the related evolutionary track on the HR diagram. On the 
other hand, if the mass loss is turned on later on the RGB phase of the progenitor star, when 
the mass of the helium core is larger (AMj^ e ~ 1%)) the final model with 0.33 Mq succeeds 
to ignite He-burning and eventually evolves to the final WD stage with a C/O core of about 
the 53% of the total mass. Figure [8] shows its evolutionary track. This is the lowest C/O WD 
we managed to produce for this chemical composition, in fact the models with mass smaller 
than 0.33 Mq do not ignite the 3a and finish their evolution as He-core WDs. The present 
evolutionary computations prove that it is possible to have C/O WDs with mass as low as 0.33 
Mq, significantly lower than 0.5 Mq, the classical and commonly accepted lower limit. This 
means that in the mass range 0.33 - 0.5 Mq both He and C/O core WDs can exist. Thus it 
is interesting to compare the main characteristics of their structure and evolution. A detailed 
description of the cooling evolution of low mass WDs has been recently published by Panei et al. 
(2007). Figure [9] shows the comparison between the core chemical profiles of the two remnants of 
0.33 Mq. While in one case the core is constituted by nearly pure helium, in the other a mixture 
of carbon and oxygen is present. Notice that in the center there is a significant abundance of 
carbon, Xi2c= 0.285, as in standard C/O WDs, thus these very low mass WDs should not be 
called simply oxygen WDs. As previously stated, a peculiar feature of these very low mass C/O 
WDs is that the helium-rich buffer is very thick, up to about 50% of the total mass, while the 
value of typical C/O WDs is of the order of 1-2 %. The thickness of the hydrogen rich outermost 
layer is almost the same in the two models, that is M= 0.0014 Mt^d for the WD with a helium 
rich core and M= 0.0015 ^iwD for that with C/O one. 
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Figure 9. Abundances in mass of 
helium, carbon and oxygen as a function 
of the mass coordinate for the two WDs 
with M= 0.33 M , the He-core (dashed 
line) and the C/O core (solid lines). 
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Figure 10. Comparison between the 
cooling curves, logL/L/0 vs. time, of the 
two WDs of 0.33 M , that with the He- 
core (dashed line) and with C/O core 
(solid line). 
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Figure 11. Radius vs. effective 
temperature of the two WDs of 0.33 M , 
that with the He-core (dashed line) and 
with C/O core (solid line). 



_ 1 1 1 1 1 1 1 

- — CO WD 
: — - He WD 


, | , , , , | , , , , j , , , , j , , , 




~ / / 
: / / 
/ / 
/ / 
= / / 

Ih / / 

: / / 
: / / 
- / / 


WD: M=0.33 M o 




: / / 
- / / 


Progenitor: M MS = 2.3 M o Z=0.04 Y=0.32 




: / / 
J / 
- / 

y 






5 4.25 


4 3.75 3.5 3.25 
log T [K] 





Figure 12. Surface gravity vs. effective 
temperature of the two WDs of 0.33 M Q , 
that with the He-core (dashed line) and 
with C/O core (solid line). 



Figure [TU] shows the comparison between the cooling curves of the He-core WD (dashed line) 
and the C/O core (solid line). As expected, the remnant with the helium rich core cools slower 
than that with a C/O one, as the specific heat per gram of helium is larger than that of carbon 
and oxygen; thus the thermal content of the He WD is larger than that of a C/O one with 
the same total mass. Such a difference in the cooling times is slightely counterbalanced at 
IoqL/Lq ~ -4 as the C/O core crystallizes, with the concomitant energy release, while the He 
core does not, in the presently computed range of temperatures and densities. 

Figures [TT1 and [T2l show the comparison between the radii and surface gravities, respectively, 
as a function of the effective temperature for the two WDs of 0.33 Mq. For a given effective 
temperature, the He-core WD is more expanded than the other one, difference not negligible 
at the beginning of the cooling (~ 8.5%). Figure [TBI shows the comparison between the central 
density p c as a function of the luminosity for the two WDs of 0.33 Mq, with the He-core (dashed 
line) and C/O one (solid line). We did not show the plot with the central temperatures T c 
versus luminosity relations for the two WDs because they are nearly identical. 
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Figure 13. Central density p c vs. 
luminosity of the two WDs of 0.33 M , 
that with the He-core (dashed line) and 
with C/O core (solid line). 



4. Conclusion 

In conclusion, as already anticipated by Iben & Tutukov (1985), we proved, by means of fully 
and consistent evolutionary computations, that the minimum mass for a C/O WD is about 
0.33 M , near the minimum possible mass of the helium-core at the onset of He-burning, which 
occurs at the RGB phase transition (M~ 2.3 M Q ). 

As a consequence, in the mass range 0.33-0.5 M Q both He and C/O core WDs can exist. As 
expected and already shown by Panei et al. (2007), the cooling times of these two classes of 
WDs are quite different, being the He-core remnants significantly slower than the C/O ones. In 
addition, we showed also that the He WDs are more expanded that their C/O counterparts at 
a given effective temperature. 
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